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We outline a generic, flexible, modular, yet efficient framework to the computation of en- 
ergies and states for general nanoscopic systems with a focus on semiconductor quantum 
dots. The approach utilizes the configuration interaction method, in principal obtaining 
all many-body correlations in the system. The approach exploits the powerful abstract- 
ing mechanisms of CH — h's template facility to produce efficient yet general code. The 
primary utility of the present approach is not in the resulting raw computational speed, 
but rather in minimizing the time from initial idea to final results. 
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1. Introduction 

The configuration interaction (CI) approach to electronic structure calculations has 
the benefit of being conceptually simple and the capability of being potentially 
exact — exact in the sense that all eigenstates and energies of the given model system 
can be exactly computed given sufficient computational resources. The detriment 
of the CI method is the great resources it requires. Generally speaking, state of 
the art calculations perform at 8-10 particles. The reason is that the CI method 
is essentially a direct diagonalization of the system's Hamiltonian matrix, and the 
size of this matrix grows (at least) exponentially with the number of particles. 

Because of this, other techniques, such as quantum Monte Carlo (QMC) and 
density functional theory, in all their various guises, are often employed (and often 
required) in place of CI calculations. Of these, QMC and CI methods are the two 
primary tools in cases where correlations are strong, where excited states in addition 
to ground states are required, and especially where the many-body states themselves 
(as opposed to the energies) are required. QMC has great potential for being far 
more efficient than CI, yet is also far more complex both in concept and in practical 
aspects of coding. Thus, QMC may be viewed as being well suited for problems 
where a large number of particles is required, or for stable production code which 
will be applied to a specific fixed set of systems over a prolonged period of time. 
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One may distinguish, however, between fast code written slowly — that is, taking 
a long time to code, test, debug, and maintain — and slow code written quickly. QMC 
belongs in the former category, while CI belongs in the latter. Particularly for basic 
research, and perhaps less so for applied research, a more useful interpretation of 
"efficiency" than the actual computation time is to include development time as 
well. Often, the most relevant time is that between initial idea and final result, and 
by no means is this duration necessarily dominated by raw computation time. In 
this sense, a properly structured CI base can be an important piece of a research 
framework, particularly when flexibility is highly desired. 

2. The need for a flexible computing environment 

Flexibility is particularly required for computations involving the basic physics of 
synthetic nanostructurcs such as semiconductor quantum dots."'^ There are four 
main reasons for this: the nonuniversality of synthetic structures, the insufhciency 
of ab initio methods for large systems, the many different environmental couplings in 
heterogeneous solid-state devices, and the highly tunable symmetries experimentally 
attainable. In this section, we explore these issues by comparing them to atomic 
systems, where very robust and highly efficient code is available. 

In semiconductor quantum dots — for example, artificial atomJ2 defined by elec- 
trostatic (lateral) confinement of electrons in a two-dimensional electron gas — each 
manufactured device is unique and, for example, the confining potential of two de- 
vices with ostensibly identical gate geometries and identical numbers of confined 
particles, will nevertheless have unique confining potentials, tunneling barriers, and 
other parameters significantly affecting the transport, electronic, and spin proper- 
ties of the device. Contrast this with atoms, where every distinct atomic species is 
an indistinguishable particle; a large-scale, robust, sophisticated, and efhcient code 
basJS eventually pays off since the resulting application is indefinitely applicable to 
all future problems of the same atomic species. 

In addition to the non-uniqueness of synthetic structures, many are also large 
in size, and can contain, for example, several millions of different atomic species 
arranged in heterogeneous geometries. Not only are current ab initio methods in- 
sufficient to deal with nonperiodic systems of such size, but often they are also 
undesirable; numerical results on simpler effective Hamiltonians many times yield 
more insightful, more intuitive results, particularly when studying experimentally 
observed effects (decoherence, for example) whose underlying physical mechanisms 
are unknown. In such cases, numerical results are more properly viewed as aids to 
theoretical analysis rather than pure simulation. 

In quantum nanoelectronic systems, and mesoscopic systems more generally, the 
relevant length scale of the constituent particles — for example, their de Broglie wave- 
length, or their mean free path — extends over the dimensions of the confinement 
potential. In such cases, boundary effects or other couplings to various environ- 
ments often drive the observed physics. These couplings generally destroy quantum 
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coherence and lead the system to more classical-like behavior. Experiments can 
now probe this boundary between classical and quantum physics and the relevant 
physics can often times be explored only numerically. A numerical approach must 
be flexible enough to investigate these quantum statistical systems, neither periodic 
nor isolated, coupled to a heterogeneous and fluctuating environment, and driven 
far from equilibrium. Such a purely quantum numerical framework does not cur- 
rently exist. Indeed even the basic physical theory is often poorly understood. In 
such cases, particularly when only a few interacting electron are present (but infi- 
nite bath degrees of freedom), flexibility usually (not always) trumps raw processing 
power. 

3. The choice of programming language 

In terms of raw speed, simpler languages are usually superior. This, along with 
simple inertia, likely explains the dominance, or at least the preponderance, of 
FORTRAN for numerical code. But, generally speaking, simpler implies less expressive 
which in turn implies less flexible. If flexibility is to be put at a premium, then a 
more expressive language is required. This will generally come at a cost of slower 
speed, but quicker development. 

A practical compromise is a language where the programmer can decide how 
much of an abstraction penalty to pay. This pay-only-for-what-you-use approach is 
well represented in the C-I--I- programming language, and this is the language in 
which we have developed the Diagqi0 framework we describe below. 

More important than efficiency, more important than flexibility, correctness is 
the single most important requirement of all numerical code. An aid, but by no 
means a savior, to correctness is C++'s type safety. A constant can really be a 
constant without recourse to preprocessor macros. Pointers can point to (various) 
constant objects or can be constant pointers themselves. Parameters to functions 
can be passed by value, by reference, or, most efficiently, by constant reference. 
These often obviate the need for raw pointers to memory and all their error prone 
complexity. Authors of a class have at their disposal (and their discretion) a great 
deal of control over how memory is precisely allocated, and how objects are created, 
assigned, copied, or converted to other types. These language features are great tools 
in producing flexible yet robust numerical frameworks. 

3.1. Generic programming 

None of the above has anything to do with object-oriented programming — the 
paradigm most frequently associated with C++. Indeed C++ supports several pro- 
gramming styles, including procedural and functional, as well as object-oriented. 
However, in our opinion, the single most useful feature of C++ in regards to nu- 
merical computation is its generic programming facilities. These can be seen as a 
bridge between expressiveness on the one hand, and efficiency on the other. Particu- 
larly for scientific programming, recent development in generic programming have 
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clearly shown that flexibility and abstraction need not necessarily incur a run-time 
penalty and have moved the state of the art far beyond simple parametrized types. 

At its most basic level, generic programming separates data types from 
algorithms;^ Thus, for example, a single sort algorithm can be written that can sort 
objects of arbitrary type, including user-defined classes, so long as the expression 
A < B IS defined for objects of those types. And because of operator overloading, 
the Icss-than operator can be defined by the author of the data type. Because of 
this, generic programming can be said to make code forwards compatible in time, 
rather than simply backwards compatible with previous versions. If a specific type 
requires a more efficient algorithm, that particular type can be made into a special 
case through partial template specialization. 

Importantly, this great fiexibility can come with little or no performance 
penalty — a crucial difference compared to dynamically typed languages such as 
Python or Ruby. The C++ generic facilities produce statically typed compiled code. 
In effect, a generic template function is a program that writes programs and is essen- 
tially a type-safe meta-programming facility. Upon compilation, the compiler takes 
a template function, and instantiates a specific version for whatever data type is 
required. Thus, a single sort template function can produce several custom-made 
(by the compiler) versions in the executable; one each for, say, integers, doubles, 
complex numbers, or many-body state vectors sorted by energy. The important fact 
is that this happens at compile time, not run time, and so the produced code can be 
very efficient. The compiler can turn this single function into an unbounded number 
of functions operating on types which the original author could not have conceived. 

The Diagon framework we describe below really only scratches the surface of 
what generic programming can do. For example, templates can be used to keep track 
of dimensional quantities,'^so that, for example, the system knows that length/time 
is a velocity and can give a (compile-time) error when a length is assigned to a 
velocity, even though all quantities are, say doubles. In fact, scientific programmers 
recognized quite early orplthe benefits of generic programming. 

As an example of the flexibility of generic programming coupled with function 
overloading, we give here an example of their use in the Diagon framework which 
we discuss more fully below. Specifically, we can consider the inner product of var- 
ious state vectors describing electrons in a quantum dot. There are at least three 
types of states that are generally required. At the single-particle level, there are 
the base single-particle orbital states |a), where a denotes a full set of quantum 
numbers, which we label generically as SPState. A many-body Fermion system 
is described by antisymmetrised vectors \ao,ai,a2, ■ ■ ■) (Slater determinants, la- 
beled AntiSymmState) in accordance with the Pauli exclusion principle. Thirdly, 
correlated states — ^^Ci\a}^,a\,a\, . . .) (induced by Coulomb interactions or 
spin symmetry ) require a description in terms of coherent superpositions of Slater 
determinants, LinCombState. 

The inner product of SPState's clearly depends on the the particular system 
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under consideration. However, the algorithmic computation of inner products of 
either LinCombState's or AntiSymmState's is essentially identical and independent 
of the underlying SPState. Thus, once supplied with an inner product function for 
the particular SPState, the calculation for inner products involving more complex 
state vectors can proceed automatically, without a need for rewriting the functions. 

The object-oriented solution to this problem is to define a class hierarchy and 
to pass pointers (or references) to the functions. At run time, the appropriate inner 
product will be called. Besides the extra (run-time) cost involved in dereferenc- 
ing objects for dynamic polymorphism, this approach will almost certainly lead to 
an ever-growing hierarchy of states, with commensurate costs in maintenance and 
testing. In addition, the core code base requires altering with every new SPState 
introduced to the system; as the number of states known to the system grows, 
maintenance, and continued testing for correctness, becomes more and more of a 
burden. 

In a generic approach, the code base remains small. Here, one defines generic 
containers AntiSymmState<SPState> and LinCombState<State>^ which can hold 
any type of state. Provided the user implements a set of SPState 's and defines 
single-particle inner-product functions, then inner products involving many-body 
states need not be (re)written; the compiler will write a custom version of the 
appropriate inner product function. The template functions themselves indicate the 
algorithm, not the data types. 

Because generic components can be combined as required, a small set of generic 
components are capable of producing combinatorially many functions. Thus, with 
a small set of components, a large array of composite objects can be defined, max- 
imizing the flexibility of the numerical approach. This flexibility need not incur a 
run-time penalty; all template instantiations are implemented at compile-time. This 
is a tremendous advantage both in the speed of developing a custom application 
from a set of generic components, as well as in maintaining efficient code at run 
time. 

We have implemented such a generic framework, Diagon, for computing eigen- 
states and eigenvalues of semiconductor quantum dots for arbitrary potentials, with 
arbitrary numbers of particles, and arbitrary Hamiltonians. This is a framework, 
meant for developing applications rather than a tool for end users. We describe 
the framework below and show that it can be used to build flexible and extensible 
applications without undue sacrifice on run-time efficiency. 

4. The Diagon Framework 

The Diagon^ framework consists of a set of generic components useful in the manip- 
ulation of many-body quantum states and computations involving them. There are 



^Note that a LinCombState can describe linear combinations of AntiSymmState's with different 
sizes. This would be required when par ticle number is not conserved, as, for example, when studying 
transport through a quantum dotll^I 
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currently components for dealing with various types of Fermion states, components 
for computing matrix elements of one and two-body operators, generators return- 
ing proper spin eigenstates given a set of singly and doubly occupied single-particle 
orbitals, and components for computing eigenstates and eigenvalues of Hermitian 
operators. These components are all generic and the single-particle states themselves 
need to be provided. We give in this section an overview of the Diagdn components. 
In the following section, we provide an example calculation for calculating spectra 
and states of two-dimensional parabolic dots with spin-orbit interactions. 



4.1. Generic quantum many-body states 

There are three classes of many-body states: AntiSymmState<SPState> is the pri- 
mary class for antisymmetrised product states. In real space, these are Slater de- 
terminants. Upon instantiation, the (generic) parameter SPState must name an 
existing class encapsulating a known single-particle state. For example, in two di- 
mensional parabolically-confined quantum dots, the single particle states may be 
the well-known Fock-Darwirffll states \mns). Typically, the inner product between 
these two single-particle states will also be provided. In the Fock-Darwin case, we 
simply have {n'm's'\nms) = Sss'^nn'^mm' ■ Operations are provided for creating and 
destroying particles in these states, as well as computing inner products if supplied 
with the inner product for the underlying single-particle states. 

In addition to the (user provided) single-particle states and the (Diagon pro- 
vided) AntiSymmState's, a generic class LinCombState<State , Coef f > is provided 
for encoding linear superpositions of states, |LinCombState) = X^a'^al'^)- Here, the 
template parameter State is the type of component states (which can be, for exam- 
ple, SPState's or AntiSymmState's) and Coef f is a template parameter encoding 
the type of the coefficients, which will usually be real or complex numbers, but, 
since the class is generic, could be of more exotic type. Operations are provided for 
adding and removing terms from the sum, for checking and setting normalization 
of the state as a whole, for indexing a particular term, and for iterating over all 
terms. If class State defines an inner product, then (LinCombState|LinCombState) 
and (LinCombState|State) are both defined. 

The final class of many-body states is the StateSet<State , Coef f > class. This 
class encapsulates a set of LinCombState<State , Coeff>'s which all share a com- 
mon basis. It is used, for example, as a return type from a diagonalization routine 
where each LinCombState is an eigenstate. It can also be used to describe a set 
of spin eigenstates. (See Sec. 14.31 1 Operations are provided, for example, to add or 
remove a basis vector to the StateSet, or to add or remove a LinCombState. 

With the above three generic components, one need only define a particular 
single-particle state and an inner product, and many-body states and linear super- 
positions of them are made available, and a full suite of inner products and other 
manipulations and computations are made available. 
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4.2. Generic operator functions 

A general quantum operator O can be given a matrix form with elements {ipi\0\iljj), 
with i,j running over all basis vectors which can be SPState's, AntiSymmState's, 
or LinCombState's. The Diagon framework provides generic operators for both one 
and two body operators. 

A general one-body operator can be written in second quantized form as 

d = Y,o,,clc, (1) 

where cj, creates a particle in state \p). In the Diagon framework, matrix elements 
of the operator O are implemented generically as 

oneBodyOpCbra, ket , matelem) . 

Here, bra and ket can be any (combination) of the three basic types of states in 
Sec. 14.11 and matelem is a user-defined function evaluating Opq in Eq. That is, 
matelem must have signature 

ReturnType matelem(SPState bra, SPState ket) , 

where ReturnType is an arbitrary type. The generic function oneBodyDp returns 
whatever matelem returns. 

We see that once the user defines single-particle properties, which will be dif- 
ferent from system to system, the Diagon framework implements the many-body 
functionality, which is in a sense a universal function of the single-particle physics. 

Two-body operators are implemented generically in a similar way. A difference 
now is how spin is treated. In particular, any spin-independent two-body operator 
{e.g., the Coulomb interaction) can be written as 

U ^ U.,jMc\„c]^,ci„,Cko, (2) 

where the Uijki = (ij\U\kl) are the coeflicient which must be provided by the 
user, with all indices i,j,k,l denoting single-particle states. This is again imple- 
mented as 

twoBodyOpCbra, ket, matelem) 

but the two states must now contain at least two particles (AntiSymmState's or 
LinCombState's) and matelem must have signature 

ReturnType matelem (SPState bral, SPState bra2, 
SPState ketl, SPState ket2) . 

The generic function twoBodyOp returns whatever matelem returns. 
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In Sec. IS] we describe a specific example using the above generic components. 
Before we do, however, we discuss generic facilities Diagon provides for computing 
eigenstates of total spin and other more general Hermitian operators. 

4.3. Generic spin states 

Because the CI technique is so computationally intensive, it is important to take 
advantage of every significant symmetry in the system as this affords a possibil- 
ity to block-diagonalize the Hamiltonian matrix, drastically reducing the computa- 
tional load. Simple symmetries such as conservation of spin or angular-momentum 
projection along a given axis {S^^^ or L^^^, say) are simple to implement since 
these symmetries do not produce correlations and their conservation can always 
be encoded in a single Slater determinant. Other symmetries — total spin, S'^qj, be- 
ing the most prominent — do induce correlations, and a single Slater determinant 
(AntiSjrmmState) cannot in general be written down in which S'tQ^ is a good quan- 
tum number. 

In such correlated basis may be used which preserves the many-body 

symmetry. This, in general, would require a prediagonalization step. However, 
the SU(2) symmetry of spin, along with its higher-dimensional representations, al- 
lows all eigenstates of spin to be written down for arbitrary orbital configuration, 
essentially relying on the appropriate products of Clebsch-Gordon coefficients. 

Such a facility is provided in Diagon through the spinGen generic function, 

spinGen(AntiSyinmState<SPState> config, int twoS, int twoSz) , 

which returns a StateSet. (See Sec. 14.11 ^ Each element of the StateSet is a 
LinCombState<AntiSyinmState<SPState>> which is an eigenstate of S'^^^ with the 
appropriate spin quantum number. Input to spinGen is two times the spin S and 
two times the projection Sz- (This is to keep the inputs integers.) Input is also the 
orbital configuration as an AntiSymmState<SPState>. 

4.4. Generic diagonalization 

The CI method eventually requires a diagonalization. Currently, Diagon employs 
the uBLAS linear algebra library of the Boost*^ project, along with a bindings 
library allowing C-I--I- to directly interface with the LAPACK algorithms. This can 
be extended to other diagonalization routines without much trouble. The function 
has signature 

vector<double> diagon (Matrix H, 

StateSet<State , Coeff> eigenVecs, 
size_t numEigs) , 



''Boost jhttp:/ /www. boost. org} provides a collection of free peer-reviewed C-l — h libraries with an 
emphasis on generics and portability. 
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where H is a uBLAS matrix. The final two arguments are optional. If the first 
is provided, then the eigenvectors of H are calculated and placed in eigenVecs. 
Otherwise, only eigenvalues are computed. The final argument numEigs indicates 
how many eigenvectors and eigenvalues to compute. If omitted, all are computed. 
The function itself returns the eigenvalues of H in a vector<double>. 

To construct H, one would normally call oneBodyOp and/or twoBodyOp for each 
of the elements. To aid in this, a function matrixDp is provided 

Matrix matrixDpCbasis , matelem) 

which returns the uBLAS matrix obtained by applying matelem fScc. 14.1^1) to each 
of the basis vectors in basis. This function calls oneBodyOp and twoBodyOp as 
appropriate. 

As mentioned above, the primary purpose of the Diagon framework is to pro- 
vide flexible, generic tools to aid development. This aspect was placed at a higher 
premium than pure computational efficiency, although the generic nature of the 
framework is very well suited to producing efficient code as well. In the following 
penultimate section, we provide an example illuminating the strengths of the Diagon 
framework. 

5. Example: Parabolic dots vifith spin-orbit interactions 

Using the generic framework Diagon, a complete diagonalization program can be 
set up and run remarkably quickly. Once the basic (non-generic) components are 
provided, the package, at compile-time, produces a custom-made set of classes and 
functions dealing with linear superpositions of many-body states. These can be used 
as a simple diagonalization to obtain spectra, or the eigenstates can be used for fur- 
ther computations in, for example, problems of quantum dynamics and decoherence, 
where the actual states are required, and where correlations in the system play an 
important role. In this section, we outline a diagonalization procedure to illustrate 
the use of the Diagon framework. 

We look specifically with at a two-dimensional quantum dot in a GaAs/AlGaAs 
heterojunction parabolically confined in the plane and in the presence of both spin- 
orbit interactions and a magnetic field perpendicular to the plane of the dot .-H" The 
Hamiltonian may be written as 

H = Hqd+Hso, (3) 

where the quantum dot Hamiltonian is given by two harmonic oscillators plus a 
Zeeman term 

ffqd - nn+ (^a^a f^^- (b^b + ^) + Qf^BB'^S', (4) 

with n± = [{cu^ + 4wg)i/2 ± l]/2. Here, ujq is the confinement frequency char- 
acterizing the parabolic confinement \^ and ujc = eB'^ / {mc) is the cyclotron fre- 
quency. The final term in Eq. is the spin-orbit interaction. We take a lin- 
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int mainC) { 
// 1. Raw parabolic parameters tCiaAs) 
double omgM = 1; // meV 
double bField = 1; // Tesla 

FDParam fiiParam - genFOPuramGaAsCoinegciS, bField); 

// Z. Raw spin-orbit parameters (GaAs) 
double rashbaAlpha = 0.1; // meV nm 
double dresselBeta = 0.6; // meV nm 

FDParamSO fdParamSO = genFDParamSOGaAsCfdPoram, rashboAlpha, dresselBeta); 

// 3. Basis vectors 

siie_t nuiiBasisStates = 1000; 

GenFDState fdGenCfdParam) ; 

vector<FDState> basisStatesCnumBasisStates) ; 

generateCbasisStates . beginC) , basisStates.endC), fdGen); 

4. Hamiltonian 
FDMatElemDiag fdMatElemCfclParam); 
Matrix homilFD = matrixOpCbasisStates, fdMatElem); 
FOMatEl emSO f dMatEl emSOCf dParamSQ) ; 

MatrixComplex htmilSO = matrixOpCbosisStates , f dMatElemSO) ; 
MatrixComplejs hanil - hamilFD + hamilSO; 

// 5. Diagonalization 

StateSet<FDState, Complex> eigenVteettrs; 

eigenvectors . addBasisCbasisStoteis.b^tnO , basisStates . endC)> ; 
vectQr<double> eigenVolues = dtagonChamil , eigenvectors, numToKeep); 

// Print the results ... 

return 0; 

} 

Fig. 1. Example cod showing minimal mainO function for computing spin-orbit eigenvalues 
and eigenvectors. There are only 21 lines of code. 



earized model including both Dresselhaus and Rashba termsf^ respectively given 
by Hso = l3(—UxPx + <^yPy) + Oi((jyPx — <JxPy), whcrc (Tfc are the Pauli matrices, and 
(3 and a are respectively the Dresselhaus and Rashba coefficients. In terms of the 
Bose operators of Eq. we can write 

Hso = AS'+ [Vl^iah^ + i[3b) - VL+{aa + i^a^)] + h.c, (5) 

where h.c. is the Hermitian conjugate, is the spin raising operator, and — 

The objective in this example is to diagonalize Eq. Q in the basis of the eigen- 
states of Eq. given by the well-known-'^-'^ Fock-Darwin states \nms). 

Figure^shows a minimal function which does this using the Diagon framework. 
The example is meant only for illustrative purposes; header files and additional 
comments have been stripped for brevity. Section 1 of the code simply creates 
parameters for Eq. Q); the genFDParamGaAs () function takes loq and the external 
field and computes VL± and the Zeeman energy in Eq. Q using GaAs material 
parameters. Section 2 performs a similar task for Eq. 

Section 3 generates the basis states with which to perform the diagonalization. 
The object fdGen () is a function which is instantiated with the fdParam object; 
successive calls of fdCenO return the Fock-Darwin state \nms) with the next- 
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highest energy. Thus caUing fdGenO 1000 times will yield the 1000 lowest-energy 
states with the given material and model parameters. The line 

generate(basisStates.begin() , basisStates . endO , f dGen) ; 

does exactly this. The basis states are then stored in the vector basisStates. 

Section 4 of Fig.n]creates the Hamiltonian matrix, Eq. Q, using the matrixDp 
function of Sec. 14.41 The function objects f dMatElem and f dMatElemSO are each 
to be called with two arguments (Fock-Darwin states, FDState) and return matrix 
elements of TJqd and Hgo respectively. Thus, these are the matelem objects described 
in Sec. EJ 

Finally, in section 5 of Fig.Q] the diagonalization occurs. A StateSet fScc. I4.l|l 
is first created and the basis states are added to it. Then, diagon (Sec. I4.4|l is 
called and the eigenvalues are placed in the vector eigenvalues, whereas all the 
eigenstates, each an orthogonal superposition of the basis vectors, are placed in the 
StateSet eigenvectors. These results can then be printed, or otherwise processed. 

We stress that this simple example is meant for illustrative purposes only. Many 
of the functions in Fig. ^ take optional arguments and support different interfaces 
and much of the Diagon framework has not been explicitly mentioned in this exam- 
ple. Nevertheless, it does illustrate how a properly constructed generic framework 
can support a flexible computing environment. For example, elliptic dots could eas- 
ily be added to the above example. One would simply need to define an additional 
matelem function object which computes the appropriate matrix elements, and an 
additional parameter object containing the eccentricities and so forth. This addi- 
tional term could then be added to the Hamiltonian through the matrixOp function. 
Additionally adding Coulomb interactions among the particles would proceed along 
essentially identical lines. 

6. Conclusions 

We have outlined an approach to the CI method which utilizes the generic pro- 
gramming facilities provided by the CH — h programming language. The general idea 
is that much of the CI machinery is independent of the actual single-particle states 
used. A generic approach allows one explicitly separate algorithms and data types 
and allows a great deal of code reuse. This has been implemented in the Diagon 
framework, which focuses on (but is by no means restricted to) two-dimensional 
semiconductor quantum dots. A generic approach such as this offers a good com- 
promise between rapid development and flexibility on the one hand, and efficient 
code on the other. 
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